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FOREWORD 


'  The  free  electron  laser  instability  for  a  solid  relativistic  electron  beam 
propagating  in  combined  transverse  helical  wiggler  and  uniform  axial  guide 
fields  is  investigated  within  the  framework  of  the  linearized  Vlasov-Maxwell 

equations.  Ht  is  assumed  that  v/yb  <<1 ,  where  v  is  Budker's  parameter  and 

2  -'v 

Ybmc  is  the  electron  energy.  'Stability  properties  are  investigated  for  the 
choice  of  equal ibri urn  distribution  function  in  which  all  electrons  have  the 
same  value  of  the  linear  combination  of  transverse  and  helical  invariants./ 

C1  -  2ybmu)bCh  =  const.,  and  a  Lorentzian  distribution  in  the  axial  invariant  Cz> 
(Here  wb  is  a  constant.)  The  instability  growth  rate  is  calculated  including  a 
determination  of  the  optimum  value  of  the  ratio  of  beam  radius  to  conducting 
wall  radius  (Rq/Rc)  for  maximum  growth.  It  is  found  that  the  maximum  growth 
rate  for  a  solid  electron  beam  is  comparable  to  that  for  a  hollow  beam  with 
similar  parameters.  Moreover,  the  introduction  of  a  small  axial  momentum  spread 
(^/Ybirar^ff  few  percent)  significantly  reduces  the  instability  growth  rate,  pr- 
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1.  INTRODUCTION 

1-9 

In  recent  years,  the  free  electron  laser  instability  has  been 
extensively  Investigated  with  particular  emphasis  on  applications 
to  intense  microwave  generation.  For  the  most  part,  previous  theoretical 

analyses  of  this  instability  have  been  carried  out  either  for  uniform 
3-6 

density  beams  with  infinite  transverse  dimension,  or  for  annular 
7  8 

electron  beams.  ’  The  present  paper  examines  the  influence  of  finite 
radial  geometry  on  the  free  electron  laser  instability  for  a  solid 
electron  beam  propagating  in  combined  transverse  helical  wlggler  and 
uniform  axial  guide  fields.  The  analysis  is  carried  out  within  the 
framework  of  the  linearized  Vlasov-Haxvell  equations,  including  a 
determination  of  the  optimum  value  of  beam  radius  R^  for  maximum 
growth  rate. 

The  present  analysis  is  carried  out  for  an  infinitely  long  relativ¬ 
istic  electron  beam  propagating  in  the  combined  transverse  wiggler  and 
uniform  axial  guide  fields  described  by  Eq.  (1).  Equilibrium  and  stability 
properties  are  calculated  for  the  specific  choice  of  electron  distribution 
function  (Eq.  (5)], 

fb  '  i(Ci  "  2Ybra"bCh  “  2Yb'nT1)G(Cz)  , 

where  n^,  y^.  and  Tx  are  positive  constants,  Ci(  C^,  and  Cz  are 

the  transverse,  helical,  and  axial  invariants*^  defined  in  Eqs. 

(6)  -  (8),  and  the  axial  distribution  function  is  normalized  according 
to  dCzG(Cz)  ■  1.  Equilibrium  properties  are  investigated  in  Sec.  II, 

and  stability  properties  are  examined  in  Secs.  Ill  and  IV,  assuming 

2 

that  v/y^  «  1,  where  v  is  Budker's  parameter,  y^mc  is  the  characteristic 

electron  energy,  and  m  is  the  electron  rest  mass.  j 
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In  Sec.  Ill,  making  use  of  the  linearized  Vlasov- Maxwell  equations, 
we  obtain  the  coupled  eigenvalue  equations  (33),  (37),  (40),  (45),  and 
(46)  that  describe  free  electron  laser  stability  properties  in 
circumstances  where  the  perturbed  transverse  fields  can  be  approximated 
by  the  vacuum  waveguide  fields.  For  short  wavelength  perturbations, 
the  axial  component  of  the  perturbed  longitudinal  field  can  be  approximated 
by  (Eq.  (49)], 


£X<'> ' 


\  0  ,  ot 


r/R0)  , 


0  <  r  <  R„ 


otherwise 


where  Jj(x)  is  the  Bessel  function  of  the  first  kind  of  order  l, 

6.  ,  is  the  s'th  zero  of  J  (B,  ,)  -  0  and  $  ,  is  a  constant.  In 

Sec.  IV,  substituting  Eq.  (49)  into  the  coupled  eigenvalue  equations, 
we  obtain  closed  algebraic  dispersion  relations  for  the  transverse 
electric  (TE)  and  transverse  magnetic  (TM)  polarizations. 


Introducing  the  normalized  dimensionless  function  [Eq.  (59)], 

2  J£(x) 

Gts-(x)  ■  S,.-  tt~ r— 2  * 

(X  6l, s’> 

it  is  shown  in  Sec.  IV  that  the  coupling  between  the  longitudinal  and 


transverse  perturbations  Is  proportional  to  G  , (o  . ,  En/R  ) 

£ S  fc' i|8  U  C 

for  the  TE  mode,  and  to  G.  Rn/R  )  for  the  TM  mode.  Here 

£s  £+l,s  0  c 

“t+l.S  and  e£+l,s  are  the  sthzeroes  °f  j;+lS+l,s)  ■  °  and  Ji+l(B£+l,s)  “  °* 
respectively,  Rc  is  the  radius  of  the  outer  conducting  wall,  and  the  prime 

denotes  (d/dx)J  (x).  Assuming  that  the  maximum  of  the  function  C  ,  (x) 

£  £S 

occurs  at  i  ■  x  , ,  we  note  that  the  maximum  instability  growth  rate 
Is 

occurs  at  a  value  of  R„/R  given  by  Rn/R  »  x.  for  the  TE 

0  c  1  0  c  is  l+l,s 

mode,  and  by  Rq/Rc  ”  x^g,/8t+^  s  Bor  t*le  ®  mode.  This  result  is 

g 

different  from  that  obtained  for  a  hollow  electron  beam. 


8 


NSWC  TR  81-145 


A  detailed  numerical  analysis  of  the  TE  mode  [Eq.  (62)]  and  TM 

mode  [Eq.  (63)]  dispersion  relations  is  presented  in  Sec.  IV.  Two 

features  are  noteworthy  from  the  numerical  analysis.  First,  for  the  optimized 

value  of  Rn/R  ,  the  instability  growth  rate  for  the  TM  mode  is  comparable 
u  c 

to  that  for  the  TE  mode.  Moreover,  the  growth  rate  is  reduced  substantially 
by  introducing  a  small  amount  of  axial  momentum  spread  (A/y^mc  £  0.01). 


9/10 


NSWC  TR  81-145 


II.  EQUILIBRIUM  THEORY  AND  BASIC  ASSUMPTIONS 


The  equilibrium  configuration  consists  of  a  relativistic  electron 
beam  propagating  in  the  combined  transverse  wiggler  and  uniform  axial 
guide  fields  described  by 


B  =  -SBcos (0  -  k0x)£r  +  6Bsin(8  -  k0z)ee  +  B^ 


where  Bq  and  6B  are  constants,  and  kg  is  the  axial  wavenumber  of  the 
helical  wiggler  field.  In  Eq.  (1),  cylindrical  polar  coordinates  (r,6,z) 
are  used,  with  z-axis  along  the  propagation  direction,  and  and 

are  unit  vectors  in  the  r-,  8-,  and  z-directions ,  respectively.  In 
the  present  analysis,  we  assume  that  the  axial  wavenumber  of  the  helical 
wiggler  field  is  sufficiently  large  that 

"l - -  i-  ~\  «  U0H  ,  (2) 

lw0  “  V  V  00 

where  wg  "  k^V^,  “c  =  e®0^Ybmc  *S  t*ie  e^ectron  cyclotron  frequency, 

2 

Rq  is  the  characteristic  beam  radius,  Y^mc  is  the  characteristic 

2  1/2 

electron  energy,  c  is  the  speed  of  light  in  vacuo,  »  c(yfa  -  1)  / y ^ 

is  the  mean  axial  velocity  of  the  electron  beam,  and  -e  and  m  are  the 
electron  charge  and  rest  mass,  respectively. 

It  is  also  assumed  that 


v/yv  «  1. 


where  v  *>  N,e  /me  is  Budker's  parameter, 


,2*  ,R_ 

d6  i 


dr  r  nb(r,  6  -  kQz)  , 


is  the  number  of  electrons  per  unit  axial  length  of  the  beam,  n^Cr,  e_  k^z) 
is  the  equilibrium  electron  density,  and  Rc  is  the  radius  of  the  conducting 
wall.  The  inequality  in  Lq.  (3)  indicates  that  the  beam  is  very  tenuous, 


...  »V:— oin  —ii<ia  r&KTmKn* 
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In  the  parameter  regimes  of  practical  interest  for  free  electron 

laser  applications,  the  axial  distribution  function  G(C^)  is  strongly 

peaked  about  «  YbmVb  "  const,»  with  characteristic  half-width 

AC  «  y  mV  .  Moreover,  in  the  present  analysis.  It  is  also  assumed 
Z  D  D 

that  the  axial  motion  is  nonresonant  with 


.2  2  ,  2 

k0vx  *  “c  ’ 


where  vz  *  pz/ym  is  t^le  axial  velocity  of  a  typical  beam  electron. 


(10) 


We  therefore  approximate  Eq .  (8)  by 


in  6B 


10 


Cz  “  Pz  +  (l  -JTbI  [prcos(e  -  ko*>  ■  P6sin(e  -  koz)]- 

0  c  0 


(11) 


Making  use  of  Eqs.  (6),  (7),  and  (11),  it  is  straightforward 


to  show  that  the  combinat  ii>o  CA  •  In  Kq.  (5)  ran  be  expressed  as 


10 


Ci  -  2Ybinch  «  (pr  -  )cos(e  -  k0*)J 

4 

+  (pe  -  YbmV  +  ^(4^)Sln(6  '  k°z)] 


(12) 


+  y2b*2+0ir,  0  -  kQz)  , 


where  tuQ  -  kQVb>  “c  “  eBo^Ybmc’  and  the  effective  potential  ^(r,  0  -  kyZ) 
is  defined  by 


fu c  "  “b\  r_  SB 
“c1^  -  wj  k0  B0 

.  ( 3LL%Y  4  f  V?  +  -  (C  -  YbnV, 

\“o  '  “c /  J  'Bo'  b  V\  “0  /  b  r  b  b 


^0(r’  6  "  k0z)  "  (u*b“c  ■  %^r  +  2 

k2  2 


sin(e  -  kQz) 


(13) 


As  a  simple  example,  we  consider  an  axial  distribution  function 


in  which  all  electrons  have  a  same  value  of  C  ,  i.e., 

z 
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where  the  mean  radius  Rq  is  defined  by 


y.  m  V  -  u  /  2  \  B„  /  b 
b  0  c  wn  0 


(“b“c  ~  V 

and  use  has  been  made  of  Eq.  (2).  Additional  general  equilibrium 
properties  associated  with  the  distribution  function  in  Eq.  (5), 


including  helical  distortions  of  the  beam  equilibrium  for  finite  6B/Bq, 
are  discussed  in  Ref.  10. 


14 
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III.  LINEARIZED  VLASOV-MAXWELL  EQUATIONS  FOR  A  TENUOUS  BEAM 


In  this  section,  we  make  use  of  the  linearized  Vlasov-Maxwell 
equations  to  investigate  the  free  electron  laser  stability  properties 
of  a  relativistic  solid  electron  beam  described  by  the  equilibrium 
distribution  function  in  Eq.  (5).  We  adopt  a  normal-mode  approach  in 
which  all  perturbations  are  assumed  to  vary  with  time  and  space 
according  to 

6\J/(x,t)  =  l  i^n\r)exp{l[£6  +  (k  +  nkQ)z  -  wt])  ,  (17) 

i,n 

where  lmm  >  0.  Here,  uj  is  the  complex  eigenf requency f  k  +  nk^  is  the 
axial  wavenumber,  and  l  and  n  are  integers.  Moreover,  it  is  also  assumed 
that  the  perturbations  are  close  to  resonance  with 


|w  -  (k  +  nkQ)Vb|  lo^,  wc  ,  (18) 

where  ujg  »  kQVb  and  u>c  «  eBQ/Ybmc. 

The  Maxwell  equations  for  the  perturbed  electric  and  magnetic  field 
amplitudes  can  be  expressed  as 


Z  »  £(x>  “  i(w/c)g(x)  , 

V  *  B(x)  -  (4i./c)J(jc)  -  i  (w/c)E(x)  , 


(19) 


where 


”  'e 

is  the  perturbed  current  density. 

r° 

fb(x,^)  *  e  dx  exp(-iut) 
J  —  oo 


}  d3P  Z  £b(?*«> 


In  Eq.  (20), 


E(x’)  + 

%  a 


v'  x  B(x') 


(20) 


(21) 
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is  the  perturbed  distribution  function,  i  -  t'  -  t,  and  the  particle 
trajectories  ^'(t1)  and  jj'(t')  satisfy  djjc'/dt'  -  and  djg'/dt'  «  -e^'  *  ^/c, 
with  "initial"  conditions  jg'(t'  ■  t)  ■  and  v‘(t'  ■  t)  » 

Within  the  context  of  Eqs.  (3)  and  (18),  the  perturbed  distribution 
function  in  Eq.  (21)  can  be  approximated  by 


fb(^t,jg)  ~  ~  ~  J  dr  exp(-iur)|2  ^ymiui(v'  •  Jg) 

-PzU’-bl))  ^2  fS  +  (*'  *  h  i)\r,  f°}  • 


(22) 


2  22  1/2 

where  y  »  (1  +  p  /m  c  )  ,  and  use  has  been  made  of  Eq.  (19).  To  lowest 


order,  the  axial  motion  of  an  electron  is  free-streaming  with 


10 


ym 


t)  . 


(23) 


Moreover,  within  the  context  of  Eq.  (18),  on  the  right-hand  side  of 
Eq.  (22)  we  retain  contributions  to  and  in  the  orbit  integral  of 
the  form8 


U'c  j  g 

v'  *  V  - —  cos(6  -  knz  -  k_.v  t )  , 

r  z  cu_  -  a)  Bn  0  0  z  ’ 

0  c  U 


(24) 


and 


z  w0  “  “c  ®0 


C  sin (6  -  kQ*  -  kovzT^ 


(25) 


Finally,  since  the  oscillatory  modulation  of  the  radial  and  azimuthal 
orbits  is  small  amplitude  (Eq.  (2)],  we  approximate 


=  r  ,  e'  =  e  , 


(26) 


in  the  arguments  of  the  perturbation  amplitudes  on  the  right-hand  side  of 
Eq.  (22). 
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Substituting  Eqs.  (23)  -  (26)  into  Eq.  (22),  we  obtain  the  perturbed 
distribution  function 


fb(^’?)  “  E  fb")e*P  UU6  +  (k  +  nkQ)z) ) 


i.n 


iec  l  eXpll't6t  (k  +  nk°)2l7X  C  £(n) 
“  u  '  <k  nk0)vz  \  “  2  2,1 


(r)  + 


e6B 


2ymc2kg  “0 


(27) 

wnere  the  function  Xn,(j>,ii,,k)  is  defined  by 


3fh  8fb 

Xn'(£"“»k)  “  2(ym4i  -  (k  +  n'kjj)pz]  — j  +  (k  +  n'kg)  - - 


(28) 


3  Pi 


and  8  *  v  /c.  In  Eq.  (27),  the  term  proportional  to  X  is  the 

z  z  n 


longitudinal  portion  of  the  perturbed  distribution  function. 

Similarly,  the  terns  proportional  to  X^+1  and  X  j  in  Eq.  (27)  are  the 
transverse  electromagnetic  portions  of  the  perturbed  distribution  function. 


Consistent  with  Eq.  (18),  the  eigenf requency  ui  can  be  approximated 


7  2  2 

by  u  =  (k  +  nko^Vb'  We  therefore  approximate  u  /c  -  (k  +  nkg  +  kg)  by 


ui2/c2  -  (k  +  nkQ  +  kg)2  » 


(k  +  nk_)  _ 

- ~  +  2k0(k  +  nkQ)  +  k2 


<  0  , 
(29) 


for  k  +  nkg  >  0.  Evidently,  Eq.  (29)  indicates  that  the  n  +  1  mode  in 
Eq.  (27)  is  a  non-propagating  wave  in  a  vacuum  waveguide.  Without  loss 
of  generality,  for  a  tenuous  beam,  we  therefore  assume 


£(n+l) 

r,i-l 


(r) 


p(n+l) 

8,1-1 


(r) 


0  . 


(30) 


in  the  subsequent  analysis.  Making  use  of  Eq.  (30),  ^(x,^)  in  Eq.  (27) 
can  then  be  expressed  as 
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iec  exp{i[t8  +  (k  +  nk^z]) 

u  ^  u  -  (k  +  nk.)v 

1  ,n  0  2 


VB  E*nJ<r) 

|  n  z  z  i  x 


+  AAn-l^r,l+l(r)  '  i^0,£+l(r)  1 


}• 


(31) 


(32) 


where  the  dimensionless  parameter  A  is  defined  by 

e6B  U0 

A  "  ,  2,  u0  -  u  ’ 

2y^mc  kQ  0  c 

and  use  has  been  made  of  the  approximation  y  -  y^,  which  is  consistent 
with  Eq.  (18). 

From  Poisson's  equation,  V  •  £(jjO  “  4rp(x),  and  the 
Maxwell  equation  (19),  we  obtain  the  differential  equation, 


[V'  +  4  -  < k  +  nk0)2]i£> 


4,i(k+nk0) 

(r)  =  - - -  (r) 


(33) 


for  the  axial  (longitudinal)  component  of  the  perturbed  elecLric 


field  £<">.  In  Eq.  (33),  pf°(r) 


-e  i 


d^pf^"^  is  the  perturbed  charge 


density,  V2  i  r-1(3/3r)  (r3/3  r)  -  i2/r2,  and  use  has  been  made  of  J^(r) 
Vbp£n^(r).  In  the  tenuous  beam  limit  (Eq.  (3)],  the  transverse  field 
components  (r)  in  Eq.  (31)  can  be  approximated  by  the  vacuum 

g 

waveguide  fields.  In  this  context,  the  present  stability  analysis 
utilises  the  vacuum  transverse  electric  (TE)  and  transverse  magnetic 
(XM)  waveguide  modes  as  a  convenient  basis  to  represent  the  general 
electromagnetic  field  perturbation  ^"^^(r),  which  is  determined  from 


^  -  (k  +  nkQ 


k0> 


(34) 


“  l1  32  fe2?£+l(^)  1  C  Iz  *  ?i®2,*+l(?) 


Making  use  of  Eqs.  (19)  and  (34),  and  neglecting  the  perturbed 
current  density,  the  vacuum  waveguide  fields  can  be  expressed  as 
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l£>>  - 


for  the  TE  mode,  and 
i(n-l) 


loR 


(35) 


b  -  J.  K+i  >  . 


“t+l.s  *+l.*Va»+l..r/Kc 


2,£+l 

Un-l) 


(r)  “  £W,sJt+l('WlSr/Rc>  • 


(3b) 


F  — -  «(«-!)..  _  .  (k  +  nk0  -  k0>Rc 

r.*+l(  6,t+l(r)  1  ^ -  Lt+l,sWl,sr/Rc>  - 

for  the  TO  mode.  In  Eqs.  (35)  and  (36),  b£+1>s  and  c£+ls  are  constants. 

Jj£ *  (x)  is  the  Bessel  function  of  first  kind  of  order  fc '  and  a 

’  fc+l,s 

and  6£+1>s  are  the  sth  roots  of  j;+1(a£+1>s)  -  0  and  3^(8^)  -  0, 
respectively.  Here  the  prime  (')  denotes  J£lJ(x)  «  (d/dx)  J£+£  (x) . 

Attei  auue  straightiorwai.1  algebraic  uianipu lat ion  of  Eqs.  (19),  (35),  and 
(36),  it  can  be  shown  that 


f  2  a2 

^ - q  i  .i  _  -,2  t+l,s 

|c2  (k+nk0'V 

i  R 

c 


b  j  ( -iM »  s  \ 
£+l,sJ£+i\  r  ) 


-£{h  !>>;"(„}  . 


(37) 


for  the  TE  mode,  and 


^  4  -  (k  +  nk.  -  k  ,)2  -  g^1»£  r  J  (  bl+1>sT  ) 
c2  0  V  r2  fX+l,sJ£+lV  R  ) 


(38) 


’  M(k  +  nko  *  ko)^+il)(r)  -  7  3^m(r))  ■ 

for  the  TO  mode.  Moreover,  making  use  of  the  continuity  equation, 

i(k  +  nkn  -  k„)j(n_1|  =  i  L.  [rj(n-l)  l  +  Ml  ±  1)  j(n-l) 
1+1  0  k0'jz,£+1  r  Jr  rJr,i+l'  +  r  Je,l+I* 

(39) 
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the  approximation  =  V^^j^fr)  (consistent  with  Eq .  (3)], 

and  approximating  k  +  nk^  =  kQ/(l  -  V^/c)  on  the  right-hand  side  of 
Eq.  (38),  we  find  that  tq .  (38)  can  be  expressed  as 


^2  -  (k  +  nkQ  - 
c 


k0) 


"  i+1 , s 
2 


£  j  /fi+iial 

t+i.s  i+n  r 


An  Is _ 

rc  |3r 


1  rJ 


(n-1) 
r,  i+1 


(r»  +  i(^)J^(r)}  . 


(40) 


for  the  TM  mode. 

For  convenience  of  notation  in  the  subsequent  analysis,  we  intro¬ 
duce  the  effective  susceptibility, 


Xn  n> (w.k)  =  4ne2  d3p 


A  ^ , u) , k) 

U)  —  (k  +  nkQ)vz 


(41) 


Moreover,  to  simplify  the  present  analysis,  we  also  assume  that  the 
beam  rotation  is  slow  with 


“b  <<  V  “o  •  (i2) 

Within  the  context  of  Eq.  (42),  we  can  show  from  Eq.  (12)  that  the 
equilibrium  distribution  function  is  an  even  function  of 


pr  -  2ybmcAcos(8  -  kgZ)  (43) 

and 

pQ  +  2YbmcAsin(6  -  kQz)  ,  (44) 


for  the  beam  rotations  satisfying  <<  wc>  ujq.  Making  use  of 
Eqs.  (31),  (41),  and  (43),  the  perturbed  charge  and  current  densities 
are  given  by 


‘(n-1) 

e,i+i 


(r) 


4irw 


Gj(u,k,r) 


- (n-1) .  . 
r,i+lU' 


icAp^n^ (r) 


(45) 


where  the  function  Gg(u>, k,r)  is  defined  by 
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C^w.k.r) 


*  6  E(n) 

*n ,  n  z  2,1 


(r) 


+  A*n,n-l[i^(r) 


^(n-1)  .  .  . 

XLe,f+llr' J 


(46) 


Equations  (33),  (37),  and  (40),  when  combined  with  Eq.  (45),  constitute 
one  of  the  principal  results  of  this  paper  and  can  be  used  to  investigate 
stability  properties  for  a  broad  range  of  system  parameters.  Moreover, 
in  limiting  cases,  the  dispersion  relation  for  the  free  electron  laser 
instability  can  be  obtained  in  a  closed  form  (Sec.  IV). 
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XV.  FREE  ELECTRON  LASER  STABILITY  PROPERTIES 


In  this  section,  simplified  expressions  are  obtained  for  the  longi¬ 
tudinal  perturbations  in  Eq.  (33), and  the  results  are  used  to  derive  the 
dispersion  relation  for  several  values  of  azimuthal  harmonic  number  t. 

The  present  analysis  assumes  short  wavelength  perturbations  with 

|q*|  =  |(k  +  nk0)2  -  u2/c2|  »  1/R2  .  (47) 

Moreover,  for  u  =  (k  +  nk^)V^  and  k  +  nk^  -  k^/(l  -  V^/c),  the  inequality 
in  Eq.  (47)  can  be  expressed  in  the  equivalent  form, 

(1  +  Vb/c)2Y^k2R2  »  1  ,  (48) 


which  is  readily  satisfied  in  the  parameter  regimes  of  present  experimental 
interest.  As  shown  in  Appendix  A,  for  short  wavelength  perturbations 
satisfying  Eq.  (48),  the  axial  component  of  the  perturbed  electric  field 
E^n ^ (r)  in  Eq.  (33)  can  be  approximated  by 

Z  ,  l 


p(n) 

z,£ 


(r) 


0  * 


e/Rq)  ,  o  _<  r 
otherwise  . 


<  K, 


(49) 


In  Eq.  (49),  6  ,  is  the  s'th  root  of  J  (G.  ■)  =  0,  and  it  ,  is  a 

St  f  S  A  1  J  S  JL  y  S 

constant. 

Substituting  Eqs.  (35)  and  (49)  into  Eq .  (46),  multiplying  Eqs. 

(33)  and  (37)  by  rJ£(6£  g,r/R0)  and  rJi+l^u£+l  s r/F'c>»  respectively, 
and  integrating  from  r  *  0  to  r  »  R^ ,  we  obtain  two  homogeneous  equations 
relating  the  perturbation  amplitudes  ,  and  b  .  For  the  TE 

S  AT  2.  y  S 

mode  polarization,  these  are 

Vf  C'  '  ‘"o  -  *>  (7  A,.  *\*^)  j2.  (\’-  ) 

b  0  (50) 
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and 


(k  +  nk_)R 
_  0  c 

Ybat+l,s 


i+l 


fR 

c 

*SJ0 


dr  r 


~T  (k  +  nk-  -  kj2  -  +  A2) 


n,n-l 


,2  Ial+  l,sr 


1+1 


c 

(51) 


irHsji q*v.-|0c  '•<«»  -  ~s^I)Jt+i(-!5^£')  • 


where  •(x)  is  the  Heaviside  step  function  defined  by 


«(x) 


/  1  •  *  ‘  0  • 

[  0  ,  otherwise 


(52) 


Similarly,  for  the  TM  mode  polarization,  we  obtain 

.  2 


and 


rR 

•*!; 


(53) 


♦.,.■1;^ «  *<»o  -  «>  (t  «...  *  t  ) 

^  »»  fc. 1 1  .  «  1  it r  . 

y*’  (  >  &  \ 

£+l\  R  ) 

c 

(54) 

B0  s'  -  fKC  /  S,r  \ 

(k  +  nk0)R0  A*t,s’j0  dr  r  «<R0  -  r)xn>nJl+1  (  f-  ) 


'  (k  +  nkQ  k0)^  - 


--  +  A  x 


n,n-l 


*Vl(  P. 


El+l,sr 


)• 


where  usa  has  been  made  of  w  =  (k  +  nk_)V,  and  (k  +  nk  )  =  k,./(l  -  V, /c). 

Ob  o  0  b 

In  the  present  analysis,  it  is  assumed  that  the  axial  distribution 
function  has  the  form 


G(C  )  =  * 


<C_  -  y b"Vb)  +  6 


2  ’ 


(55) 


where  A  is  the  characteristic  spread  in  C  about  the  mean  value  C  «  y,  mV,  . 

*  z  b  b 
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We  further  assume  that  the  characteristic  spread  A  is  small  in  comparison 
with  YbmVb.  Substituting  Eqs.  (5)  and  (55)  into  Eqs.  (28)  and  (41), 
we  obtain  the  approximate  expression 


“2  "  +  nk0Hk  +  n’k0)c2 


n,n 


Yb^o  -  U  +  nkQ)Vb  +  i|k  +  nk0|i/Ybm]2 
0  .  Rq  <  r  _<  Rc 


0  <  r  '  V 
(56) 


In  obtaining  Eq.  (56),  use  has  been  made  of  Lc .  (18).  Making  use  of  the 

definition  of  Budker's  parameter  in  Eqs.  (3)  and  (4),  the  term  4v/y,R2 

b  0 

n  2  2  2 

in  Eq.  (56)  can  also  be  expressed  as  4v/y,  R“  «  to  ,/c  ,  where  w  ,  « 

b  U  pb  pb 

2 

4nnge  /v^m  is  the  plasma  frequency-squared. 

Tlie  condition  for  a  nontrivial  solution  to  Eqs.  (50)  and  (51) 

is  that  the  determinant  of  the  coefficients  c„  .  and  b  be  equal 

t,s  i+l,s 

to  zero.  After  some  algebraic  manipulation,  we  find  that  the  TF.  mode 
dispersion  relation  can  be  expressed  as 


ui  -  (k  +  nkg)Vb  +  i 


|k  +  nk0|a 


{4-  <k  *  "»o  -  V2  - 


If  |k  +  nk„  |  A 

2  \ 

(  “  -  (k  *  nk0)Vb  +  1  3 

-  4  VC  ) 
3_2  / 

V  Yba 

YbR0  J 

4A 


q  A 

2  vc 


YrR2 
b  c 


kQ(k  +  nkQ  -  kjj) 


VU 


P. 


-,E  lat+l,sR0 
££ss'\ 


^tss-f^r1-0)!"-  <k+nk0)Vb+i 


|k  +  nkQ | A 


3 

Ybm 


-4^H  '( 

3_2  Iss  \ 

YbR0 


a£-tl  ,SR0 


where  the  coupling  coefficient  QR  ,  (a..,  Rn/K  )  is  defined  by 

K.SS  i|S  v  C 


*1+1. s 


'ft+l.s 


(£  +  1)2  j2+lS+l,s> 


and  the  functions  Gist(*)  and  Hffls,(x)  are  defined  by 


(57) 


(58) 
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2  J?(x) 

C£s'(x)  ’  7.J—3 


and 


(x" 

,2 


6£.s‘> 


ll£ss'(x) 


Wx)  -  Jt(x)Jt^2(x) 
Gts,(x) 


(59) 


(60) 


In  Eq.  (57),  the  subscript  s  and  s'  represent  the  radial  mode  numbers 
of  the  transverse  and  longitudinal  perturbations,  respectively. 

For  small  viggler  amplitude  (\  «  1),  we  investigate  free  electron 

laser  stability  properties  for  u  and  k  +  nkg  near  the  simultaneous 

2  2  2 

zeros  of  the  transverse  dispersion  relation,  u.  -  (k  +  nkp  -  k^)  c  - 
2  2  2 

ai!.+l  sC  /Rc  =  anc*  t*le  ^on8iCu^ina1  dispersion  relation 


w  -  (k  +  nkp)V^  +  i 


|k  +  nkQ|i 


y.m 


-  4 


v3R2 

YbR0 


0  . 


(61) 


In  this  regard,  making  use  of  Eq.  (61),  the  TE  mode  dispersion  relation 
in  Eq.  (57)  can  be  approximated  by 


(  H—  -  (k  +  nk  -  k  ) 2  -  ■■  ~ a  1 
|2  v  0  0;  2  I 

l  c  R  J 

c 


k  +  nk 


it)  —  (k  +  nk^)V^  +  i 


o^l ; 

) 


-  7. 


3  2 
YbR0 


44 


2  vc‘ 


Y.R 

b  c 


kQ(k  +  nk0  -  kQ)  -  2 


0E  /  °H-1 ,  sKQ 
^iss’V  R 


In  a  similar  manner,  from  Eqs.  (53)  and  (54),  we  obtain  the 


c 

(62) 


approximate  TM  mode  dispersion  relation, 

„2  1 


{: 


2 

2  ” 

-  4 


8 


(k  +  nkp  -  kQ) 


£+l)S\ 


-  (k  +  nkQ)Vb  +  1 


|k  +  nkQ | 4 


1 2 


Ybm 


2  1 


3  2 

Vo 


4A 


2  vc 


kQ(k  +  nkjj  -  kp) 


"t+l.s 


M  / 
Q£ss’  ( 


(63) 
Bi+l,sKQ 


where  the  TK  mode  coupling  coefficient  Q,  , (B„ Rn/R  )  Is  defined  by 
r  0  £ss  £+l,s  0  c 

(64) 


<ss'(x)  "  G£s'(x)/J£+2{6£+1,s>  ' 


and  the  function  Gfs,(x)  is  defined  in  Eq .  (59). 
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Figure  1  shows  plots  of  C^s,(x)  versus  x  obtained  from  Eq.  (59) 

for  (a)  s’  “  ®o  i  and  (*>)  gt  “  Except  in  the  case 

6.  i  =  gn  ,  ,  tlie  plots  of  C  ,  (x)  for  arbitrary  6  ,  are  similar 

to  those  for  g^  s'  ”  ^1  3  *n  ^*6'  ^(h).  As  shown  in  Fig.  1(b),  the 

quantities  cf  .  £  G.  , (x.  ,)  and  x.  ,  denote  the  maximum  value  of 
Is  is  is  is 

Gts,(*)  and  the  corresponding  value  of  x  for  a  specified  g^  ^ , .  For 
example,  in  Fig.  1,  (x^,,  G*gl)  =  (0,  0.09)  for  BfL  g,  »  gQ  1  and 
(x0  ,  >  G?  i )  “  (9.8,  0.064)  for  g  ,  »  g  Shown  in  Fig.  2  are 

X-s  £S  £  ,  S  1  ,  J 

plots  of  (a)  X£.  1  an<^  0>)  tlie  corresponding  values  of  g|  ,  = 

G  , (xis,  )  for  several  different  values  of  the  azimuthal  and 

radial  mode  numbers  i  and  s'.  It  is  evident  from  Fig.  2(b)  that 

C^s,  decreases  rapidly  with  increasing  values  of  the  mode  numbers  i 

and  s'.  Moreover,  we  note  from  Fig.  2(a)  that  x{c_,  can  be  approximated  b 


Xis'  =  6i.s-  '  S’  y  1 


(65) 


In  this  regard,  for  s'  y  1,  C*  ,  can  he  approximated  by 

K.S 

CL'  -  Ws-*  -  Ji+i(fci,s->- s’  ' l-  (66) 

r 

Shown  in  Fig.  3  are  plots  of  (a)  Q  , /G  ,  for  the  TE  mode 

£SS  £K 

and  (b)  Q^ss,/G^s,  for  the  TM  mode,  obtained  from  Eqs.  (58)  and  (64) 
respectively.  Note  that  the  curves  in  Fig.  3  are  independent  of  tiie 
longitudinal  radial  mode  number  s'.  Evidently,  tlie  values  of  Q  ,/G  , 

fSS  £S 

and  Q  ,/G  ,  increase  with  increasing  values  of  azimuthal  and 

£ss  £s 

transverse  radial  mode  numbers,  l  and  s.  After  careful  examination 
of  Eqs.  (58)  and  (64),  we  find  that  the  maximum  coupling  between  the 


transverse  and  longitudinal  modes  occurs  for  a  value  of  R^/K^  given  by 
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is' 

1 

t+l,s 

,  /g 

is' 

1  pi+l,s 

R0/Rc  '  \  "  *■**■  (6?) 

Equation  (67)  is  valid  only  when  x^,  _<  g  for  the  TE  mode,  and 

xis-  -<fW  for  the  ™  raode-  For  xis'  >  Vi.s  (TE>-  or  xt8-  >  6i+l.B 
(TM),  the  maximum  coupling  occurs  for  RQ/Rc  “  1.  The  maximum  coupling 

coefficients  corresponding  to  Eq.  (67)  can  be  determined  from  Figs.  2(b) 

and  3.  For  example,  for  (i,s,s')  -  (3,2,1),  we  determine  that  the 

maximum  coupling  coefficient  and  the  corresponding  ratio  Rq/R  •  are 

given  by  (QEgs,  rq/rc)  “  (1.607,  0.625)  for  the  TE  mode,  and  (Q*'ggl,  Ro/Rc) 

(1.83,  0.52)  for  the  TM  mode. 

It  is  instructive  to  examine  the  present  results  for  perturbations 

with  the  lowest  mode  numbers,  i.e.,  (i,s,s')  «  (0,1,1),  particularly 

for  a  beam-filled  waveguide  with  Rq/Rc.  ~  l.  In  this  limit,  from  Fig.  1(a), 

we  obtain  G„  , (a  )  *■  0.4  for  the  TE  mode,  and  G„  ,(B,  ,)  “  0.045  mode. 

Is  1,1  is  1,1 

We  therefore  conclude  that  the  TE  mode  polarisation  is  the  most  unstable. 
Multiplying  G  , (u.  ,)  -  0.4  by  QE, /G,  -  4.2  in  Fig.  3(a),  the 

£S  1  «  X  JLSS  fcS 

E 

coupling  coefficient  is  given  by  Qq^  “  1.7.  Assuming  zero  axial 

2  2 

momentum  spread  (a  *=  0)  and  short  axial  wavelengths  (k^R^  >:> 
the  TE  mode  dispersion  relation  in  Eq.  (62)  can  be  approximated  by 


.  2  \ 

\  *  \ 

I  ^2  +  nk0  ’  ] 

<[w  “  (k  +  nkp)V^]2 - 3 

V  c  ' 

l  *b  1 

,2  2  .2 
ulpbk0  ’ 


(68) 


for  the  U.s.s’)  =  (0,1,1)  perturbation  and  R./R 

0  c 


1.  Equation  (68) 


is  similar  in  form  to  the  result  obtained  by  Davidson  and  Uhm  for 


a  uniform  density  beam  with  infinite  cross  section.  In  particular, 
the  constant  numerical  factor  on  the  right-hand  side  of  Eq.  (68)  is 
equal  to  3.4,  whereas  in  Ref.  3  the  constant  numerical  factor  is  equal  to  8. 
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Finally,  we  have  investigated  detailed  stability  properties 
by  solving  the  dispersion  relations  in  Eqs.  (62)  and  (63)  numerically 
for  a  broad  range  of  system  parameters.  Defining  the  normalized 
Doppler-shifted  eigenf requency  by 

ii  -  [w  -  (k  +  nk0)Vfc)/k0c  ,  (69) 

we  calculate  the  normalized  growth  rate  0^  =  lmO  from  Eqs.  (62)  and 

(63).  Shown  in  Fig.  4  are  plots  of  the  normalized  growth  rate 

versus  (k  +  nkQ)/k0  for  U,s,s')  =  (3,2,1),  k^  =  10,  yb  *  10,  v/yb  “ 

2 

0.02,  and  A  •  0.01,  with  (a)  R../R  =  x,,/u.  .  for  the  TE  mode,  and 

U  c  jl  4,2 

(b)  Rn/K  «  x.,/6,  -  for  the  TM  mode.  For  these  optimized  choices  of 
U  c  31  4,2 

R0^Rc’  t*1e  instabilicy  growth  rate  for  the  TM  mode  is  comparable  to 

that  for  the  TF.  mode.  Moreover,  the  growth  rate  is  reduced  substantially 

by  introducing  a  small  amount  of  axial  momentum  spread  (A/ybmc  £  0.01). 

We  conclude  this  section  by  pointing  out  two  areas  in  which  the 
analysis  can  be  extended.  First,  the  restriction  to  very  short  ;wave- 
length  perturbations  [Eq.  (46)]  can  be  removed  in  a  relatively  straight¬ 
forward  manner.  Second,  paralleling  the  self-consistent  theoretical 

g 

formalism  developed  in  previous  studies,  the  stability  analysis  can 
also  be  carried  out  without  making  the  approximation  that  the  transverse 
perturbations  are  represented  by  the  vacuum  waveguide  fields. 
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V.  CONCLESIONS 

In  this  paper,  we  have  examined  the  free  electron  laser  instability 

for  a  solid  relativistic  electron  beam  propagating  in  the  combined 

transverse  wiggler  and  uniform  axial  guide  field  given  in  Eo.  (1). 

The  analysis  was  carried  out  within  the  framework  of  the  linearized 

Vlasov-Maxwell  equations.  The  equilibrium  (Sec.  II)  and  stability  (Secs. 

Ill  and  IV)  properties  were  investigated  in  detail  for  the  choice  of 

distribution  function  in  which  all  electrons  have  the  same  value  of 

the  linear  combination  of  transverse  and  helical  invariants,  C,  -  2v,nu),  C,  , 

b  b  h 

and  a  Lorentzian  distribution  in  the  axial  invariant  [Eqs.  (5)  and 

(18)].  One  of  the  most  important  conclusions  of  this  analysis  is  that 

the  maximum  instability  growth  rate  for  a  solid  electron  beam  is  comparable 

to  that  of  i  hollow  beam  win  sin-ilai  parameters.  horeover,  it  is  also 

found  that  the  maximum  growth  rate  occurs  at  a  value  of  Rq/R^  corresponding 

to  Rn/R_  *  x  ,  /a . , ,  for  TE  mode  perturbations,  and  F.,,/R  «  x„ 

^  c  £s  t+l,s  0  c  Js  J+l,s 

for  Tfl  mode  perturbations.  For  these  optimized  values  of  R./R  ,  the 

U  c 

instability  growth  rate  for  the  TK  inode  is  comparable  to  that  for  the  TL 
mode.  Moreover,  the  growth  rate  is  substantially  reduced  by  introducing 
a  small  amount  of  axial  momentum  spread  (i/y^mc  ^  0.01). 
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FIGURE  2  PLOTS  OF  (a)  xgs-  AND  <b)  THE  CORRESPONDING  G$s>  =  Ggs-  (xgs<)  FOR  SEVERAL  VALUES 
OF  AZIMUTHAL  AND  RAIDAL  MODE  NUMBERS,  i.  AND  S’. 
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(b)  TM  mode 
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FIGURE  3  PLOTS  OF  (a)  Q|s*/Ggs-  AND  (b)  Q^s'/Ges' 

EQS.  (58)  AND  (64)  FOR  SEVERAL  VALUES  OF  8  AND  S. 
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FIGURE  4(a)  PLOTS  OF  NORMALIZED  GRWTH  RATE  VERSUS  (k  +  nk0)/ko  FOR  (e,8.s’>  -  (3,2,1),  7b  “  10, 
vlyb  -  0.02,  AND  A2  =  0.01,  WITH  (a)  Rq/R,.  -  x31/c*4  2  F0R  THE  TE  M0DE- 


FIGURE  4(b)  R0/Rc  =  x31/04  2  FOR  THE  TM  MODE 
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APPENDIX  A 

LONGITUDINAL  PERTURBATIONS  FOR  THE  FREE  ELECTRON  LASER  INSTABILITY 


In  this  Appendix,  we  investigate  properties  of  the  longitudinal 
perturbations  about  an  electron  beam  propagating  through  a  cylindrical 
waveguide  with  radius  R^.  In  the  present  analysis,  it  is  assumed  that 
the  perturbations  have  short  wavelength  with 

U*1  -  | <k  +  nkQ)2  -  w2/c2|  »  1/Rq  ,  (A.l) 

which  can  also  be  expressed  as 

qnR0  "  (1  +  V'^b^O  ~  1  *  (A‘2) 


for  the  frequencies  of  interest  for  free  electron  laser  applications. 
Equation  (A. 2)  is  easily  satisfied  in  parameter  regimes  of  present 
experimental  interest.  In  the  limit  of  a  small  uiggler  amplitude  (A  -*  0) , 
we  obtain  the  longitudinal  eigenvalue  equation. 


/  1  3  3 _  l  2  \  ‘  (n)  .  , 

(  r  3r  r  3r  “  2  ~  qn  )  Ez,i(r) 


3 _,2  ’ 


[w  -  (k  +  nk0)Vfa  +  ijk  +  nk0|d/ibm] 
from  Eqs.  (33),  (45),  (46),  and  (56).  In  Eq .  (A. 3),  C(x) 


step  function  defined  in  Eq.  (52),  and  w 


pb 


4vc2/YbKo  is 


frequency-squared. 

For  notational  simplicity,  we  define 


(A. 3) 


is  the  Heaviside 
the  plasma- 


i/(r)  =  E<"j(r)  . 

Z  |  X 


(A.  4) 


Inside  the  electron  beam  (0  <  r  <  Rq)  >  Eq-  (A.  3)  can  be  expressed  as' 


12 
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{$hrh--2  +  r2)6*t(l)  ■ 0  •  °-‘r  <Ro  • 


(A. 5) 


where 


2  .  2 
T  =  q 


2  /  2 
Mpb/Yb 


3_,2 


"  [[“  -  (k  +  nk0)Vb  +  i | k  +  nk0|A/vbm] 


1>.  (A. 6) 


Outside  the  electron  beam  (R.  <  r  <  R  ),  Eq-  (A. 3)  reduces  to  the  free- 

u  c 

space  eigenvalue  equation 

(rS7  r  IT-4  -  *'n)4*t(r)  *  °  ■  R0  ‘  r  <Rc  ;  (A‘7> 


The  solution  to  Eq.  (A. 5)  that  remains  finite  at  r  ”  0  is 


6^n(r)  -  itJt(TR)  ,  0  r  <  RQ  ,  (A.  8) 

where  J^(x)  is  the  Bessel  function  of  the  first  kind  of  order  £,  and 

*  2  2 

$  is  a  constant.  Noting  a  R'  >>  1  in  Eq .  (A.  2),  »’e  can  express  the 
£  'n  U 

solution  to  Eq.  (A. 7)  as 


6^ut(r)  -  C[l£(qnr)K£(qnRc)  -  K  (q^)  I£  (q^)  ] ,  ^  <  r  <  »C,(A.9) 

where  1^  and  are  modified  Bessel  functions  of  order  i  ,  and  C  is  a 
constant. 

The  boundary  conditions  on  iif^r)  at  the  surface  of  the  electron 
beam  are  given  by 

■  K«w0  •  <*•“> 

and 

l£)'Sr)4»‘in3r.R0  ■  *  (A’U) 

from  Eq.  (A. 3).  Substituting  Eqs.  (A. 8)  and  (A. 9)  into  Eqs.  (A. 10)  and 
(A. 11)  gives 


A-2 


NSWC  TR  81-145 


TR 


j;<tr0> 

0  J£(tr0) 


h(%) 


(A. 12) 

-  p  1i(,.W,»Rc)  '  Wc>Wo> 

'  V°  V'nW^V  ’  VqnPc)K£(qnR0>  ' 
where  the  "prime"  notation  denotes  derivative  with  respect  to  the  complete 
argument  of  the  Bessel  function,  e.g.,  J'(TR0)  -  [dj  (x)/dx)  .  The 

£  U  £  X~TRq 

expression  for  the  longitudinal  wave  admittance  h(qn)  in  Eq.  (A. 12) 

can  be  simplified  in  several  limiting  cases,  including  short  wavelength 
2  2 

perturbations  with  |<JnRQ  |  >>  1.  In  this  case,  h(q^)  can  be  approximated  by 


h(qn)  =  "qnR0COthqn(Rc  "  V  ’  (A'13) 

and  Eq.  (A. 12)  reduces  to 
J 1 (TR  ) 

-TR0  TTfR^j"  “  qnR0COthqn(Rc  "  R0)  '  (A.1A) 

Evidently,  the  right-hand  side  of  Eq.(A.14)is  a  very  large  positive 
number,  and  the  lovest-order  longitudinal  dispersion  relation  (for  t  ->  0) 
can  be  approximated  by 


VTR0)  -■ 0  ’ 


(A. 15) 


where  T  is  defined  in  Eq.  (A. 6).  It  follows  from  Eq .  (A. 15)  that 


2  2  2 

T  R0  “  6t,s’  *  s’  "  l’  2 . 


(A.  16) 


where  fi  ,  is  the  s’th  aero  of  J  (x)  «=  0.  In  this  regard,  Eqs.  (A.  8) 
£  *  S  £ 

and  (A. 9)  can  be  approximated  by 


>£,s,J£(B£,s,r/V  ’  °  -  r  *  R0  ’ 


otherwise  , 


(A. 17) 


where  $  ,  is  a  constant.  Substituting  Eq.  (A. 6)  into  Eq.  (A. 16) 

£*  8 
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and  making  use  of  Eq.  (A. 2),  we  obtain  the  longitudinal  dispersion  relation, 


2  2 


u  -  (k  +  nk„) V,  +  i 
U  b 


|lc  +  nkfili'l  w  , 

- _  _E^  .  o 

3  2 

*hn 


(A.18) 


where  the  term  proportional  to  6*  g,  has  been  neglected. 
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